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Abstract 

Background: Exposure to ozone activates innate immune function and causes neutrophilic (PMN) airway 
inflammation that in some individuals is robustly elevated. The interplay between immuno-inflammatory function 
and genomic signaling in those with heightened inflammatory responsiveness to ozone is not well understood. 

Objectives: Determine baseline predictors and post exposure discriminators for the immuno-inflammatory 
response to ozone in inflammatory responsive adult volunteers. 

Methods: Sputum induction was performed on 27 individuals before and after a two hour chamber exposure to 
0.4 ppm ozone. Subjects were classified as inflammatory responders or non-responders to ozone based on their 
PMN response. Innate immune function, inflammatory cell and cytokine modulation and transcriptional signaling 
pathways were measured in sputum. 

Results: Post exposure, responders showed activated innate immune function (CD16: 31,004 MFI vs 8988 MFI; 
CD1 1b: 44,986 MFI vs 24,770 MFI; CD80: 2236 MFI vs 1506 MFI; IL-8: 37,603 pg/ml vs 2828 pg/ml; and IL-1 (3: 
1380 pg/ml vs 318 pg/ml) with muted signaling of immune cell trafficking pathways. In contrast, non-responders 
displayed decreased innate immune activity (CD16, CD80; phagocytosis: 2 particles/PMN vs 4 particles/PMN) post 
exposure that was accompanied by a heightened signaling of immune cell trafficking pathways. 

Conclusions: Inflammatory responsive and non responsive individuals to ozone show an inverse relationship 
between immune cell trafficking and immuno-inflammatory functional responses to ozone. These distinct genomic 
signatures may further our understanding about ozone-induced morbidity in individuals with different levels of 
inflammatory responsiveness. 

Keywords: Air pollution, Environment, Ozone, Gene expression, Human sputum, Immune response, 
Innate immunity, Systems biology 



Background 

Ozone (03) is an important criteria air pollutant known 
to induce exacerbation of asthma and other airway dis- 
eases [1,2]. Human challenge studies have shown that 03 
exposure causes an immediate nociceptive decrease in 
lung function [3], increased airway hyperresponsiveness 
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[3] and enhanced response to allergens in allergic volun- 
teers [4,5]. 03 also induces airway inflammation in healthy 
and asthmatic volunteers that is marked by inflammatory 
cell (i.e. neutrophil (PMN)) influx to the airways [6,7]. In 
healthy volunteers, the inflammatory response in the air- 
ways can occur at 03 levels near (0.08 ppm) or below 
(0.06 ppm) the National Ambient Air Quality Standard 
(NAAQS) for ambient 03 with noted variability across 
individuals, such that some are very responsive while 
others have a very weak inflammatory response, suggest- 
ing the presence of a inflammatory response phenotype in 
humans [8,9]. Previous studies in healthy adults [10] as 
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well as in animals [11] have suggested that certain genes, 
such as the GSTMlnull genotype, can act as genetic 
modifiers of the inflammatory response to 03 and may 
explain some of the individual variation often observed. 
At present, the role inflammatory responsiveness plays in 
03 susceptibility is unclear, and its association with 
ozone-induced changes in innate immune function and 
host defense genomic signatures remain largely unknown. 

In this study, our goal was to characterize factors in 
the airways that may contribute to the differential re- 
sponse of adult volunteers to ozone. Various cellular fea- 
tures associated with ozone responsiveness in the 
airways were assessed including the inflammatory re- 
sponse, innate immune function and transcriptional sig- 
naling both before and after ozone exposure. We 
stratified a cohort of 13 healthy, 4 atopic and 10 allergic 
asthmatic subjects on the basis of their PMN response 
to 03 into a responder and non-responder group. The 
inflammatory response was defined as the difference in 
% PMNs pre versus post exposure, and a PMN response 
of >10% successfully distinguished responders from non- 
responders as previously reported [12-14]. Our analysis 
focused on establishing cell-associated differences in 
sputum at baseline and post exposure in inflammatory 
responders and non-responders. Cells recovered from 
sputum derive from the lumen of the large central air- 
ways and thus reflect pathophysiologic cellular events 
in the large airways [15]. Specifically, we analyzed (1) 
genomic signatures, (2) biological indicators of innate 
immune function and (3) cellular and fluid phase mar- 
kers of inflammation in order to determine whether 
these indicators were constitutively predictive at base- 
line and/or discriminatory post exposure of the inflam- 
matory responder phenotype in adult volunteers. 
Furthermore, as it is unclear whether asthma or atopy 
status is associated with the inflammatory response 
phenotype, we examined the relationship of asthma and 
atopy status in our cohort of responders and non- 
responders. 

We demonstrate that ozone responders and non- 
responders can be differentiated both at baseline and 
post exposure using various cellular indicators. Most 
strikingly, we find that ozone-induced activation of in- 
nate immune function and inflammation is negatively 
associated with genomic signaling of immune cell traf- 
ficking networks. Specifically, responders and non 
responders have opposite genomic responses in terms of 
activation of immune cell trafficking. In addition, we 
show that responders compared to non responders have 
elevated expression of pro-inflammatory cytokines IL-8 
and II- lb at baseline, suggesting a primed inflammatory 
state, but simultaneously possess muted innate immune 
function. Asthma or atopy status per se does not influ- 
ence the inflammatory response phenotype. 



Methods 

Subjects 

Twenty seven adult subjects (13 healthy, 4 allergic non 
asthmatics, 10 allergic asthmatics), aged 21-35 years, 
with no history of smoking in the past 10 years, com- 
pleted the study (Table 1). All subjects underwent a 
physical examination, a routine blood panel with 
complete blood cell count, and differential and allergy 
skin testing for common allergens (cat, cockroach, dust, 
grass, mite, shellfish, tree, weeds). Healthy subjects were 
required to have a negative methacholine challenge test 
result, and asthma status was confirmed by clinical his- 
tory and positive (>8 mg/ml) methacholine challenge 
test. No asthmatics were receiving oral or inhaled cor- 
ticosteroid therapy but all used albuterol on an as 
needed basis. Female subjects were required to have a 
negative urine pregnancy test result before challenge, 
and all volunteers were required to be free of chronic 
cardiovascular or respiratory illness (excluding asth- 
matics) and be free of acute respiratory illness (for asth- 
matics, no exacerbations requiring oral or inhaled 
corticosteroids) within 4 weeks of 03 challenge. All sub- 
jects had FEVi and forced vital capacity (FVC) values of 
80% or greater of predicted value and FEVyFVC ratios 
of 75% or greater of predicted normal value for height 
and age. All subjects were screened for their ability to 
provide an adequate induced sputum sample during 
their training sessions 48 hours prior to exposure and 
this sample acted as their pre-exposure baseline. Pre and 
post sputum samples were separated by a minimum of 
48 hours. The study protocol was approved by the 
Institutional Review Board at the University of North 
Carolina Medical School in Chapel Hill and the US 
Environmental Protection Agency (EPA), and informed 
consent was obtained from all subjects before their par- 
ticipation in the study. 

Ozone (03) exposure and testing protocol 

The 03 exposures were conducted in an exposure 
chamber at the EPA Human Studies Facility on the 
campus of the University of North Carolina (Chapel Hill, 
NC). Each subject was exposed to 0.4 ppm 03 for two 
hours while performing four 15 -minute sessions of inter- 
mittent moderate exercise (expiratory minute ventila- 
tion, 30-40 L/min) on a treadmill, separated by 
15 minutes of seated rest. Lung function, breath sounds, 
and vital signs were assessed immediately before and 
after exposure. 

Sputum collection and analysis 

Sputum was obtained 5 hours after exposure and pro- 
cessed as previously described [10]. Sputum samples 
were used as 03 preferentially absorbs in the large cen- 
tral airways during exercise [16], the same region from 
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Table 1 Subject demographic data 



Responders 


Height (cm) 


Weight (kg) 


Age (yr) 


Gender 


Race 


Smoker 


Diagnosis 


Minute ventilation (Vmin) 


01 


157 


61.2 


21 


Female 


Caucasion 


no 


healthy 


36.8 


02 


195 


81.8 


20 


Male 


Caucasion 


no 


healthy 


29.5 


03 


171 


53.8 


30 


Male 


Caucasion 


no 


healthy 


28.3 


04 


180 


83 


31 


Male 


Asian 


no 


healthy 


27.3 


05 


182 


78.6 


23 


Male 


Caucasion 


no 


healthy 


37.3 


06 


170 


58.5 


28 


Female 


Caucasion 


no 


healthy 


35.5 


07 


153 


48.7 


22 


Female 


Caucasion 


no 


healthy 


30 


08 


159 


75.4 


22 


Female 


Caucasion 


no 


healthy 


32 


09 


186 


66.6 


23 


Male 


Caucasion 


no 


atopic + ve skin test 


32.8 


10 


167 


77.4 


35 


Male 


African American 


no 


atopic + ve skin test 


34 


11 


173.5 


82.5 


22 


Male 


Caucasion 


no 


atopic + ve skin test 


31.3 


12 


168 


92.3 


39 


Female 


Caucasion 


no 


atopic asthmatic 


34.3 


13 


181 


75.1 


20 


Male 


allergic asthmatic 


no 


atopic asthmatic 


33.3 


14 


168 


97.3 


20 


Female 


allergic asthmatic 


no 


atopic asthmatic 


30.4 


15 


178 


74.2 


31 


Female 


Caucasion 


no 


atopic asthmatic 


35 


16 


167 


69.3 


27 


Female 


Caucasion 


no 


atopic asthmatic 


32.5 


17 


170 


52 


22 


Female 


Caucasion 


no 


atopic asthmatic 


30.8 


18 


165 


73.4 


25 


Female 


Caucasion 


no 


atopic asthmatic 


31.5 


Mean (+/- SEM) 


172 (3) 


72 (3) 


26 (1) 










33 (1) 


Non-Responders 


















01 


193 


94.8 


26 


Male 


Caucasion 


no 


healthy 


32.1 


02 


156 


56.3 


27 


Female 


Caucasion 


no 


healthy 


30.8 


03 


162 


59.1 


25 


Female 


Caucasion 


no 


healthy 


34.8 


04 


172.25 


70.9 


24 


Female 


Caucasion 


no 


healthy 


33.8 


05 


157 


77.9 


26 


Female 


Caucasion 


no 


healthy 


28.3 


06 


158 


68.4 


20 


Male 


Caucasion 


no 


atopic + ve skin test 


32.3 


07 


155 


60.6 


25 


Female 


Caucasion 


no 


atopic asthmatic 


33 


08 


173 


64.8 


27 


Male 


Causasion 


no 


atopic asthmatic 


31.8 


09 


178 


64 


20 


Male 


Caucasion 


no 


atopic asthmatic 


33.5 


Mean (+/- SEM) 


167 (4) 


69 (4) 


24 (1) 










32 (1) 



which sputum samples originate [15]. In brief, three 
7-minute inhalation periods of nebulized hypertonic sa- 
line (3%, 4%, and 5%; UltraNeb 99 ultrasonic nebulizer; 
DeVilbiss, Jackson, TN) were followed by expectoration 
of sputum into a sterile specimen cup. Cell-rich "plug" 
material was selected from the raw sample and treated 
with a dilute (0.1%) solution of dithiothreitol (Sputolysin; 
Calbiochem, San Diego, CA) in Dulbecco PBS. Aliquots 
of sputum supernatant were collected and immediately 
frozen and stored at -80°C for future cytokine measure- 
ment, and the cell pellet was analyzed for cell viability 
(trypan blue exclusion), total cell count and differential 
leukocyte counts (Hema-Stain-3; Fisher Scientific, 
Hampton, NH). A portion of the remaining cells were 
also used for cell surface immunophenotyping (CD lib, 



mCD14, HLA-DR, CD16, CD80, CD86) as assessed by 
flow cytometry as previously described [10]. Acquired 
sputum samples considered acceptable for processing 
had a minimum of 75 mg of selected plug material, cell 
viability greater than 50%, and squamous epithelial cells 
less than 40%. All sample processing, including flow 
cytometric assays and slide preparations, were per- 
formed on the same day of collection. Flow analysis was 
performed on fixed samples (0.5% paraformaldehyde) 
within 24 hours of acquisition on a BD LSR-II flow cyt- 
ometer (BD Immunocytometry Systems, San Jose, CA). 

Cytokines were measured using multi-plex technology 
(Meso ScaleDiscovery/MSD, Gaitherburg, MD). Each 
sample was analyzed with the Human MIP-1 alpha Ultra 
Sensitive Kit (lot no. K0031370) and the Human TH1/ 
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TH2 10-Plex Ultra Sensitive Kit (lot no. K0031431). All 
supernatant samples were diluted 1:4 and had a final 
dithiothreitol (DTT) concentration of less than 1 mmol/L 
where no deleterious effects from the DTT have been 
observed with the MSD platform. 

Functional assays for oxidative burst activity and 
phagocytosis of opsonized zymosan were also assessed 
on sputum inflammatory cells by means of flow cytome- 
try as previously described [10]. Several phagocytosis 
parameters were assessed including the percentage of 
internalized (vs externally adhered) zymosan particles 
using Trypan Blue stain as the quenching reagent. 

Classification of the ozone induced inflammatory 
response 

Subjects demonstrated a range of responsiveness with 
respect to airway neutrophil proportion (%) following 
ozone exposure (Figure 1A). Inflammatory responsive- 
ness was defined as the % PMN post -% PMN pre 03 
exposure. In general, individuals could be grouped as 
those with a minimal inflammatory response or "non- 
responders" (N = 9) and those with a robust inflamma- 
tory response or "responders" (N = 18). An inflammatory 
response (% PMN post -% PMN pre 03) > 12 was used 
to distinguish inflammatory responders and non- 
responders and is a similar criteria as used previously by 
our group and others [12-14]. The mean ozone "re- 
sponse" for all endpoints is defined as the mean of the 
sum of the individual subject differences between post 
and pre exposure. 

FEV1 response 

Spirometry was performed on a 10.2-L dry seal digital 
spirometer interfaced to a computer (SensorMedics 
Model 1022; SensorMedics; Palm Springs, CA). At least 
three sets of qualified data were obtained and the largest 
value selected for FEVi and FVC as per American Thor- 
acic Society guidelines [17]. Pulmonary function on all 
subjects was measured using one dedicated spirometer 
and by one certified pulmonary function technician to 



minimize variability. Measurements were performed be- 
fore and immediately after exposure for use in endpoint 
analysis. 

RNA Labeling 

Total sputum RNA was prepared with Trizol (Invitrogen, 
Carlsbad, CA) homogenization and QiagenRNeasy 
column purification (Valencia, CA) according to manu- 
facturer s protocols. RNA integrity was assessed by com- 
paring 18S and 28S rRNA species using a Bioanalyzer 
(Agilent Technologies, Inc., Santa Clara, CA). RNAs 
from each group were labeled by amino-allyl nucleotide 
incorporation and reverse transcription with SuperScrip- 
tlll (Invitrogen) according to manufacturers protocol 
(Epicentre Biotechnologies, Madison, WI). Probes were 
then coupled to fluorescent-labeled cy dyes (GE Health- 
care, Piscataway, NJ), such that the experimental sam- 
ples were Cy5 (red) labeled while the human 
universal reference RNA (Stratagene/Agilent) sample 
was Cy3 (green) labeled. 

Microarray analysis 

Each Cy5 sample aRNA probe was paired with Cy3 
reference aRNA probe and hybridized to microarrays in 
formamide-containing buffer overnight at 42°C using a 
Maui Hybridization System (BioMicro Systems Inc., Salt 
Lake City, UT). The custom-designed microarray con- 
sisted of 70-mer oligonucleotides representing more 
than 2000 immune response gene targets. Commercially 
available and custom-designed oligos were purchased 
from Operon Biotechnologies, Inc. (Huntsville, AL) and 
printed onto glass slides by the Duke Microarray Facility 
(Duke University, Durham, NC). Each microarray con- 
tained targets replicate-spotted four times. The gene tar- 
get pool also included housekeeping genes (i.e. ACTB, 
GAPD, TUBA, ribosomal proteins RPL13A, RPS9) used 
as internal references, negative background controls and 
known "spiked-in" RNA oligos (Applied Biosystems, 
Austin, TX). Technical replicate assays were performed. 
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Figure 1 Ozone response for neutrophils differs between responders and non-responders. (A) Inflammatory non-responders and 
responders are distinguished on the basis of their PMN response (% PMN post -% PMN pre exposure) using an ozone response cutoff level of less 
than 12 for non-responders and greater than 12 for responders. (B) The pre versus post exposure % PMN levels differ significantly (p<0.05) for 
responders with ozone exposure, as well as the responders versus non-responders post ozone exposure. 
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Microarray data analysis 

The hybridized spot intensities were acquired using a 
GenePix 4000B scanner (Molecular Devices Corp, 
Sunnyvale, CA) and assessed using GenePix Pro 6.0 soft- 
ware. Median pixel intensities (with background sub- 
tracted) were uploaded in GeneSpring 10 (Agilent 
Technologies, Inc., Santa Clara, CA), normalized by glo- 
bal LOWESS fit [18], log base 2 transformed, filtered 
so as to exceed the background and spot quality criteria 
(i.e."Absent" probes were not included) and merged into 
datasets. Gene lists and hierarchical clusters of differen- 
tially expressed genes were generated. The use of R soft- 
ware through Bioconductor was examined to analyze 
hybridization signal intensities [19,20]. To control for 
multiple comparison, we used false discovery rate (FDR), 
estimated by 5000 permutations for each gene. 

Network analysis and enriched canonical pathways 

Molecular network and enriched canonical pathway 
analysis was performed using the Ingenuity database 
(Ingenuity Systems, www.ingenuity.com, Redwood City, 
CA). The Ingenuity database provides a collection of 
gene to phenotype associations, molecular interactions, 
regulatory events, and chemical knowledge accumulated 
to develop a global molecular network. The lists of dif- 
ferentially expressed mRNAs were overlaid onto this glo- 
bal molecular network, where protein networks strongly 
associated with the targets were algorithmically con- 
structed based on connectivity. Statistical significance of 
each network was calculated using a Fisher s exact test. 
This test generated a p value representing the intercon- 
nectivity between genes, indicating the likelihood of the 
gene products interacting by chance alone. Enriched ca- 
nonical pathways within these networks were also identi- 
fied, where a Fischers exact test was used to identify 
significantly associated canonical pathways. Canonical 
pathways with p < 0.01 were assessed. 

Statistical methods 

Baseline responses and 03 responses for all phenotype 
endpoints were compared using a Students T-test (using 
Prism version 5.01). In cases where data were not nor- 
mally distributed, non-parametric tests (Mann Whitney) 
were used. The relationship between 03-induced inflam- 
matory response and asthma status was analyzed using 
an exact chi-square test and the Mantel-Haenszel chi- 
square using S.A.S. (S.A.S, version 9.2). A p value of 
<0.05 was considered statistically significant. 

Results 

Inflammatory (PMN) response phenotype: ozone (03) 
responders and non-responders 

Sputum samples were collected from 27 volunteers 
exposed to 0.4 ppm 03 for two hours. Differential cell 



count analysis showed non-responders (N = 9) to have a 
minimal PMN response of 0.6% (±3%), while responders 
(N = 18) had a significant (p = 0.0001) PMN response of 
35% (±3%) (Figure 1A). Post exposure, PMN levels were 
significantly (p = 0.03) higher in responders compared to 
non-responders as expected, while baseline % PMN 
levels did not differ between inflammatory responders 
and non-responders (Figure IB). We note, of the 27 total 
subjects, 10 were asthmatic and 4 were atopic without 
asthma. The inflammatory response however was not 
associated with either atopy or asthma as measured by 
chi square analysis (p = 0.69). This was further confirmed 
when the subjects' neutrophil response to ozone were 
stratified by disease group and showed no difference 
when analyzed by ANOVA (p = 0.9). 

FEV1 response in inflammatory responders 
and non-responders 

Percent predicted FEV1 (%) and the "% predicted FEV1 
response" were measured to determine whether lung 
function was different at baseline and post ozone expos- 
ure in inflammatory responders and non responders. 
The mean FEV1 "response" was calculated as the mean 
of the sum of the individual subject differences between 
post and pre % predicted FEV1. We report that the % 
predicted FEV1 was not significantly different between 
inflammatory responders and non responders, respect- 
ively at pre exposure (106 ± 4% vs 107 ± 4%, p > 0.05) and 
post exposure (91 ±4% vs 93 ±5%, p>0.05). The % pre- 
dicted FEV1 "response" was also not significantly differ- 
ent between responders and non responders (data not 
shown). 

Specific cytokines increased at baseline and 
post exposure 

In order to characterize cytokine responses to 03 expos- 
ure in inflammatory responders and non-responders, 
the levels of eight different cytokines were measured 
pre and post exposure (Table 2). At pre exposure, IL-8 
was significantly elevated in responders compared to 
non-responders. At post exposure, IL-1(3 and IL-8 were 
significantly increased in responders compared to non- 
responders. Responders also demonstrated a significantly 
greater increase from pre exposure in IL-8 and IL-lp 
than non-responders (data not shown). All other cyto- 
kines (IL-2, IL-6, TNFa, INFy, IL-12) remained at rela- 
tively stable levels pre or post exposure for both 
responders and non-responders. 

Modified cell surface phenotype expression in 
inflammatory responders and non-responders 

Expression, measured as mean fluorescence intensity 
(MFI), of cell surface phenotypes on sputum monocytes, 
macrophages, neutrophils and dendritic cells was assessed 
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Table 2 Cytokine levels pre and post ozone exposure 


Cytokine (Pg/ml) 


Responders pre exposure 


Responders post exposure 


Non responders pre exposure 


Non responders post exposure 


IL-ip 


989 (511) 


1380 (672)** 


510 (184) 


318 (95) 


IL-8 


20234(10347)* 


37603 (12024)** 


3161 (1567) 


2828 (1076) 


IL-6 


104 (38) 


304 (89) 


343 (152) 


348 (138) 


IL-10 


28 (14) 


100 (62) 


55 (37) 


16(9) 


TNFa 


223 (142) 


141 (49) 


44 (14) 


63 (32) 


INFy 


204 (56) 


272 (74) 


105 (37) 


138 (51) 


IL-12 


186 (67) 


210 (104) 


209 (87) 


133 (44) 


IL-2 


203 (50) 


219 (69) 


440 (164) 


305 (101) 



Mean cytokine levels are shown (+/- standard error). 
Pg/ml = picograms per mililiter. 

* p<0.05 Responder vs Non Responder pre-exposure; ** p<0.05 Responder vs Non Responder post-exposure Pre Exposure = 48 hours pre exposure; 
Post Exposure = 5 hours post exposure. 



in inflammatory responders and non-responders pre and 
post 03 exposure. The mean 03 "response" for cell sur- 
face phenotypes was calculated as the mean of the sum of 
the individual subject differences between post and pre ex- 
posure MFI values (Table 3). Pre exposure, there were sig- 
nificant (p < 0.05) differences in surface marker expression 
between responders and non-responders for CD 14 and 
CD80 on monocytes and CD86 on macrophages. Post ex- 
posure, significant differences in expression between 
responders and non responders were observed for CD 16, 
CD80, HLA-DR and CDllb (trend) on macrophages, and 
HLA-DR on monocytes. When the mean 03 "response" 
was compared between responders and non-responders, 
significant (p<0.05) differences were found for macro- 
phage CDllb expression and monocyte CD80 expression. 
No other significant differences between responders and 
non-responders were found at pre or post exposure or for 
the calculated 03 "response" for the other surface pheno- 
types measured. 

Oxidative burst activity is not changed by ozone 

Oxidative burst activity was evaluated pre and post 03 
exposure using a functional assay. 03-induced changes 
in oxidative burst activity were not significantly different 
between inflammatory responders and non-responders 
(Table 4). Furthermore, baseline (pre 03) oxidative burst 
activity did not significantly differ between responders 
and non-responders (Table 4). 



phagocytosis indices were significantly (p < 0.05) different 
between inflammatory responders and non-responders 
at baseline in both macrophages and neutrophils. Spe- 
cifically, responders had a decrease in all phagocytosis 
endpoints in macrophages and lower particle internal- 
ization in neutrophils compared to non-responders. 
Responders did however, show higher total phagocytic 
capacity (MFI) in neutrophils suggesting greater particle 
adherence to the neutrophil surface but importantly, 
a decreased capacity to engulf and internalize the 
particle. Post exposure, neutrophil phagocytosis was 
significantly increased in responders versus non 
responders for total phagocytic capacity and for the 
number of particles taken up per neutrophil (Table 5). 
The mean 03 "response" (mean of the sum of the 
individual subject differences in phagocytosis between 
post and pre exposure) by macrophages was not 
different between responders and non-responders but 
did differsignificantly for neutrophils across several 
indices. As shown in Table 5, responders' neutrophils had 
significantly increased internalization capacity following 
03 exposure compared to non-responders who showed 
decreased internalization and a decrease in the percentage 
of neutrophils undergoing phagocytosis. Altogether, 
phagocytosis activity greatly differed between inflamma- 
tory responders and non-responders at baseline for both 
macrophages and neutrophils, while neutrophils were 
specifically affected post 03. 



Phagocytosis activity varies between inflammatory 
responders and non-responders 

Phagocytosis, assessed in sputum macrophages and neu- 
trophils, was expressed as total phagocytic capacity 
(MFI), the percentage of cells undergoing phagocytosis, 
the percentage of internalized particles and the internal 
to external ratio (I:E ratio) of particle uptake by cells 
(Table 5). These indicators were evaluated pre and post 
ozone exposure. As shown in Table 5, several 



Large differences in genomic response between 
inflammatory responders and non-responders 

A total of 140 genes were identified with significantly 
altered expression following 03 (see Additional file 1). 
Of these, 120 genes were modulated in the non- 
responder group, 13 genes were modulated in the re- 
sponder group, and an additional 7 genes were changed 
in both groups (Figure 2A). The majority of the altered 
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Table 3 Cell surface phenotype expression pre and post ozone exposure 



Surface 


R pre ozone 


R post 


R ozone 


NR pre 


NR post 


NR ozone 


Pre 


Post 


Ozone 


marker 




ozone 


response* 


ozone 


ozone 


response* 


exposure 


exposure 


response 


(MFI) 














R vs NR 


R vs NR 


R vs NR 


CDIIb Mac 


23964 (4448) 


44986 (9539) 


110 (42) 


37861 (11770) 


24770 (4280) 


8 (30) 


ns 


P = 0.07 


p = 0.05 


CD14 Mono 


5776 (559) 


11800 (1082) 


149 (41) 


9886 (1016) 


19348 (4647) 


92 (36) 


p = 0.004 


ns 


ns 


HLA-DR Mono 


7399 (849) 


10611 (832) 


85 (35) 


9715 (1412) 


15361 (2110) 


88 (45) 


ns 


P = 0.02 


ns 


HLA-DR Mac 


46144 (3293) 


48527 (5232) 


9(13) 


60440 (15219) 


76424 (9254) 


56 (29) 


ns 


P = 0.02 


ns 


CD16 Mac 


21259 (6465) 


31004 (8624) 


100 (38) 


8030 (2266) 


8988 (2011) 


27 (30) 


ns 


P = 0.01 


ns 


CD80 Mono 


2956 (603) 


2236 (576) 


16 (49) 


825 (370) 


1506 (414) 


648 (572) 


p = 0.01 


ns 


p = 0.01 


CD80 Mac 


3202 (611) 


5178 (1116) 


1 75 (88) 


3440 (1642) 


2014(649) 


32 (84) 


ns 


P = 0.02 


ns 


CD86 Mar 




J 1 QQf, (1A11) 
z. i yy\j \Z-^Z-Z-j 


1 72 (1 93) 






1 3 OT\ 

I D 


p = 0.05 






CDIIb Mono 


31267 (4186) 


32517 (3789) 


34 (25) 


27031 (3957) 


35935 (6247) 


47 (25) 


ns 


ns 


ns 


CDIIb PMN 


48975 (3856) 


53935 (6220) 


18 (16) 


51808 (8185) 


59310 (9073) 


27 (25) 


ns 


ns 


ns 


CD14 Mac 


7937 (975) 


9516 (973) 


66 (50) 


18093 (10437) 


12001 (2838) 


9(14) 


ns 


ns 


ns 


HLA-DR DC 


27304 (2787) 


27230 (2087) 


6(8) 


25986 (3583) 


29416 (2101) 


28 (25) 


ns 


ns 


ns 


CD16 Mono 


9422 (1763) 


7048 (926) 


11 (21) 


9511 (4755) 


6801 (1847) 


15 (27) 


ns 


ns 


ns 


CD16 PMN 


64160 (10852) 


62617 (9331) 


23 (25) 


65491 (29381) 


81816 (26611) 


65 (30) 


ns 


ns 


ns 


CD86 Mono 


6768 (893) 


6835 (374) 


36 (25) 


6973 (895) 


10588 (10588) 


75 (61) 


ns 


ns 


ns 



Mean (+/- standard error) values are shown for each cell surface marker, along with the corresponding cell type. 

R = Responders; NR = Non-Responders; MFI = mean fluorescence intensity; Mono = monocyte; Mac = macrophage; PMN = Polymorphonuclear neutrophil; 
ns = not significant. 

* Ozone Response values were calculated as the mean of the sum of the individual subject differences in expression (MFI) between post and pre exposure, 
pre = 48 hours pre exposure; post = 5 hours post exposure. 



genes had increased expression levels after 03 exposure 
(Figure 2B). 

To gain further insight into whether the alterations in 
gene expression profiles occurred simply as a result of 
changes in cell populations, the lists of ozone-modulated 
genes in each group were compared to a previously pub- 
lished list of 411 genes that have been suggested to be 
granulocyte-specific in comparison to other types of leu- 
kocytes [21]. Comparing the 13 genes modulated by 
ozone in the inflammatory responder group to the 
granulocyte-specific genes revealed one gene, neutrophil 
cytosolic factor 2 (NCF2), which may have been mea- 
sured at increased expression levels resulting from the 
increased presence of neutrophils found in the re- 
sponder group. However, this gene also showed 
increased expression in the non-responder group which 

Table 4 Oxidative burst activity pre and post ozone exposure 



did not show significant changes in neutrophil cell 
counts. This comparison gives evidence to suggest that 
the changes in genomic signatures induced by ozone do 
not reflect changes in neutrophil cell count. 

Biological networks potentially affected by 03 expos- 
ure were constructed using the genes that showed 
altered expression levels in responders and non-respon- 
ders. Nine networks were identified using the genes 
modified by 03 exposure in non-responders, while only 
two networks were identified in the responders (see 
Additional files 2 and 3). Figure 3 A and 3B show the 
biological networks with the highest statistical signifi- 
cance for inflammatory non-responders and responders, 
respectively. While no genes are shared between the two 
most significant networks, both contain the common 
transcriptional regulator nuclear factor of kappa light 



Oxidative burst 
activity (MFI/%) 


R pre 
ozone 


R post 
ozone 


R ozone 
response 


NR pre 
ozone 


NR post 
ozone 


NR ozone 
response* 


Pre 

exposure 
R vs NR 


Post 

exposure 
R vs NR 


Ozone 
response 
R vs NR 


Mac 


118 (16) 


84(15) 


-19 (15) 


114 (25) 


128 (31) 


116 (112) 


ns 


ns 


ns 


% Mac 


55 (8) 


36 (7) 


-20 (20) 


52 (12) 


36 (13) 


-18 (29) 


ns 


ns 


ns 


PMN 


41 (7) 


35 (5) 


14 (23) 


49 (11) 


63 (13) 


127 (134) 


ns 


ns 


ns 


% PMN 


73 (8) 


58 (7) 


35 (58) 


65 (11) 


49 (17) 


-4 (37) 


ns 


ns 


ns 



Mean (+/- standard error) values are displayed for burst activity (MFI) and % of cells undergoing burst response. 

R = Responders; NR = Non-Responders; MFI = mean fluorescence intensity; Mac = macrophages; PMN = Polymorphonuclear neutrophil; ns = not significant. 

* Ozone Response values were calculated as the mean of the sum of the individual subject differences in oxyburst values between post and pre exposure, 
pre = 48 hours pre exposure; post = 5 hours post exposure. 
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Table 5 Phagocytosis activity pre and post ozone exposure 



Phagocytosis R pre R post R ozone NR pre NR post NR ozone Pre Post Ozone 

ozone ozone response* ozone ozone response* exposure exposure response 

R vs NR R vs NR R vs NR 



Mac (MFI) 


443 (140) 


660 (138) 


139 (59) 


412 (167) 


555 (234) 


91 (50) 


ns 




ns 


ns 




% Mac 


34 (4) 


35 (4) 


17(15) 


51 (6) 


44 (8) 


-14(13) 


P = 


0.03 


ns 


ns 




Particles/Mac 


3(1) 


4(1) 


138 (59) 


3(1) 


4(2) 


91 (50) 


ns 




ns 


ns 




% Internalized Mac 


31 (5) 


32 (4) 


26 (19) 


64(10) 


57(17) 


-19 (21) 


P = 


0.02 


ns 


ns 




InternahExternal 
Particle Ratio Mac 


0.6 (0.2) 


0.5 (0.1) 


40 (27) 


2.1 (0.6) 


2.3 (0.9) 


-7 (42) 


P = 


0.02 


ns 


ns 




PMN (MFI) 


565 (144) 


542 (105) 


61 (41) 


197 (62) 


226 (27) 


78 (62) 


P = 


0.03 


p = 0.01 


ns 




% PMN 


26 (4) 


40 (6) 


70 (26) 


43 (5) 


33 (2) 


-18 (14) 


P = 


0.01 


ns 


P = 


0.01 


Particles/PMN 


4(1) 


4(1) 


66 (44) 


1 (0.3) 


2 (0.2) 


93 (58) 


P = 


0.03 


p = 0.02 


ns 




% Internalized PMN 


29 (4) 


46 (5) 


95 (37) 


72(15) 


65 (14) 


-5 (13) 


P = 


0.03 


ns 


P = 


0.03 


Internal: External 
Particle Ratio PMN 


0.5 (0.1) 


1 (0.2) 


196 (73) 


6(3) 


3(2) 


-6 (45) 


P = 


0.03 


ns 


P = 


0.04 



Mean (+/- standard error) values are displayed for each phagocytic indicator. 

R = Responders; NR = Non-Responders; MFI = mean fluorescence intensity; Mac = macrophage; PMN = polymorphonuclear neutrophil; ns = not significant. 

* Ozone Response values were calculated as the mean of the sum of the individual subject differences in phagocytosis values between post and pre exposure, 
pre = 48 hours pre exposure; post = 5 hours post exposure. 




Figure 2 Ozone exposure causes a varied genomic response 
between inflammatory responders and non-responders. (A) 

Human ozone exposure induced significant changes in the 
expression of genes, largely within the NR group. (B) A heat map 
illustrates the 140 ozone-modulated genes identified as differentially 
expressed in at least one subject group exposed to ozone. Blue 
indicates negative fold change (FC) in expression (post ozone 
challenge / baseline), and red indicates positive FC. 



chain gene enhancer in B cells (NF-kB). An examination 
of the functional role of the integrated proteins identifies 
a difference in the number of genes that encode immune 
cell trafficking-associated proteins. Specifically, the non- 
responders showed increased expression levels of 
7 immune cell trafficking-associated genes while respon- 
ders had a muted response with only 1 immune cell 
tracking- associated gene (Figure 3A,B). In particular, 
ubiquitin-conjugating enzyme E2N, (UBE2N), v-rel reticu- 
loendotheliosis viral oncogene homolog B, (RELB), beta-2- 
microglobulin, (B2M), major histocompatibility complex, 
class I, C, (HLA-C), leukocyte immunoglobulin-like recep- 
tor, subfamily B (with TM and ITIM domains), member 1, 
(LILRB1), fragment of IgG, high affinity la, receptor 
(CD64), (FCGRlAFc), fragment of IgG, low affinity Illb, 
receptor (CD16b), (FCGR3BFc), and cyclin Dl, (CCND1) 
are genes that encode proteins involved in immune cell 
trafficking. It should also be noted that FCGR3B is a gene 
that is highly expressed in neutrophils [21]. Taken to- 
gether, inflammatory responders displayed activated mar- 
kers of inflammation that included IL-8, Il-lb and 
increased levels of PMN, enhanced immune cell function 
and suppressed expression of immune cell trafficking tran- 
scripts. In contrast, non-responders demonstrated un- 
changed markers of inflammation, activated immune cell 
function and up-regulated immune cell trafficking tran- 
scripts (Figure 4). 

Discussion 

We set out to identify whether cellular and biochemical 
indicators in sputum collected both at baseline and post 
exposure could predict and characterize the inter- 
individual differences in adult volunteers who were 
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Figure 3 Differences in genomic signaling networks between inflammatory responders and non-responders. The most significant 
networks associated with ozone-induced gene expression changes in (A) inflammatory non-responders and (B) inflammatory responders are 
displayed. Networks are shown with symbols representing products of genes that are up-regulated (red symbols) by ozone, down-regulated (blue 
symbols) by ozone, or are associated (clear symbols) with the differentially expressed genes. Genes that encode proteins that play a role in 
immune cell trafficking are designated with an asterisk (*). 



grouped according to their level of inflammatory respon- 
siveness to ozone. Specifically, we examined whether 
markers of innate immune function and targeted tran- 
scripts for immune cell trafficking were able to predict 
and or discriminate between individuals with a robust 
(responder) or minimal (non-responder) inflammatoryr- 
esponse to inhaled 03. 

We report that in a group of 27 adult volunteers, 18 
were robust inflammatory responders with a mean 03- 
induced increase in % PMNs of 35%, and 9 were non- 
responders with a mean 03-induced increase in % 
PMNs of 0.6%. Responders also had significantly 
elevated levels of pro-inflammatory cytokines IL-8 and 
IL-1|3. At baseline, neither % PMN levels nor gene expres- 
sion profiles differed between inflammatory responders 
and non-responders. This was somewhat surprising since 
we have shown previously that baseline gene expression 
can be used as an indicator of inter-individual responsive- 
ness to environmental contaminants [22]. One explanation 
could be the limited selection of genes that were analyzed 
on the targeted microarray. It is possible that a whole 
genome approach would uncover baseline gene predictors 
of inflammatory responsiveness to 03. 

A select number of biological indicators did however 
distinguish the inflammatory phenotypes at baseline. For 
the responders, these included significantly decreased 
expression of cell surface phenotypes (CD 14, CD86) 



decreased phagocytosis and elevated IL-8 levels, suggest- 
ing that at baseline, individuals with the inflammatory 
responsive phenotype have muted innate immune func- 
tion in concert with a primed potential toward PMN in- 
flux due to high pre-existing levels of IL-8. 

Post exposure, responders had enhanced inflammation 
as demonstrated by increased PMN, IL-8 and IL-lp 
levels, and activated immune function in concert with 
minimal immune cell trafficking-associated genomic sig- 
naling. Interestingly, non responders had the opposite 



Ozone Exposure 



I 



Non-Responders 

No change in IL8 
No change in IL1B 
No change in PMN 

Decreased innate 
immune function 

Robust activation of 
immune cell trafficking 
pathways 



1 



Responders 

Activated IL8 
Activated I LIB 
Increased PMN 

Increased innate 
immune function 

Minimal activation of 
immune cell trafficking 
pathways 



Figure 4 Summary of the ozone-induced changes in 
non-responders and responders. 
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post- exposure responses demonstrating unchanged in- 
flammation, down-regulated immune function but with 
a robust immune cell trafficking genomic response. 
These data demonstrate the presence of a negative asso- 
ciation between a specific innate immune genomic re- 
sponse and immuno-inflammatory functional responses 
in the airways. Specifically, up-regulation of genes asso- 
ciated with immune cell trafficking appears to be nega- 
tively associated with macrophage expression of CD 16/ 
FcRylll and CDllb/CR3, both of which mediate innate 
immune responses, and neutrophil phagocytosis. It is 
important to acknowledge that these measurements 
were conducted in sputum cells that reflect events in the 
large versus distal airways, a compartmentalized re- 
sponse to ozone. Inflammatory responses are noted to 
be different in large versus distal airways. 

Since our cohort of subjects included asthmatics, we 
assessed whether disease status conferred an increased 
likelihood of being an inflammatory responder or non- 
responder. We found that neither asthma nor atopy 
(without asthma) significantly predisposed one toward a 
particular inflammatory phenotype. This suggests that 
although 03 is associated with increased exacerbation 
of asthma and can induce greater 03-induced lung 
function decline in asthmatics versus non-asthmatics 
[23,24], inflammatory responsiveness to 03 may not be 
wholly dependent on asthma status or atopy, but when 
it occurs in asthmatics, may be a risk factor for asthma 
exacerbation. 

Following 03 exposure, there were some common 
genomic responses between responders and non respon- 
ders. These included genes that are associated with in- 
flammation, specifically members of the NF-kB pathway. 
These data support research that has shown that 03 is 
known to cause increases in NF-kB binding activity in 
macrophages and type II epithelial cells [25]. Interest- 
ingly, while both groups modulated these inflammation- 
associated genes, it was the 03 non-responders who 
showed an enrichment for genomic signaling related to 
immune cell trafficking. Finally, these results identify 
through a systems biology approach, novel baseline bio- 
logical indicators and post exposure discriminators of in- 
flammatory responsiveness to 03. Our data highlight the 
complex interplay between innate immune-associated 
genomic signaling and immuno-inflammatory functional 
responses to ozone, factors that likely influence disease 
susceptibility in inflammatory responsive individuals. 

It is important to acknowledge that with these intri- 
guing data, certain limitations were present in this study. 
Most noteworthy, was the mixed sputum cell population 
that was used for the gene expression data. Absent a sin- 
gle cell population from which to draw conclusions, the 
effects of other cell types contributing to the overall gen- 
etic signaling cannot be ignored, but in this study for 



example, neutrophil levels were not likely responsible for 
the observed differences in genomic signatures. 

Conclusions 

Inflammatory responsive and non responsive individuals 
to 03 each demonstrate an inverse relationship between 
innate immune genomic signaling and immuno- 
inflammatory responses to 03, but do so in opposite 
directions. Inflammatory responders have enhanced 
immuno-inflammatory function with muted innate im- 
mune genomic signaling, while non responders have 
decreased immuno-inflammatory function with robust in- 
nate immune genomic signaling. These distinctions, par- 
ticularly in genomic signaling, may help elucidate the 
mechanisms underlying 03-induced morbidity in indivi- 
duals with differing levels of inflammatory responsiveness 
to 03. 

Additional files 

f s 

Additional file 1: Genes significantly differentially expressed after 
ozone challenge. 

Additional file 2: Networks constructed using the inflammatory 
non-responders' ozone-associated genes. 

Additional file 3: Networks constructed using the inflammatory 
responders' ozone-associated genes. 
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